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We study the isotropization of a liomogeneous, strongly coupled, non-Abelian plasma by means of 
its gravity dual. We compare the time evolution of a large number of initially anisotropic states 
as determined, on the one hand, by the full non-linear Einstein's equations and, on the other, by 
the Einstein's equations linearized around the final equilibrium state. The linear approximation 
works remarkably well even for states that exhibit large anisotropics. For example, it predicts with 
a 20% accuracy the isotropization time, which is of order t^^a ^ ^/T, with T the final equilibrium 
temperature. We comment on possible extensions to less symmetric situations. 



1. Introduction. The discovery that the quark-gluon 
plasma (QGP) behaves as a strongly coupled fluid |l] has 
provided a strong motivation for understanding the dy- 
namics of strongly coupled non-Abelian plasmas through 
the gauge/gravity duality [5]. Much has been learned 
about the (near) equilibrium properties of these plasmas 
by employing their dual description as a (slightly per- 
turbed) static black hole (see i3j and references therein). 

The more challenging phase of formation and far-from- 
equilibrium evolution of a plasma corresponds on the 
gravity side to the formation of a far-from-equilibrium 
black hole horizon and its subsequent relaxation. An 
outstanding problem in this context is understanding the 
short isotropization time of the QGP. In principle, study- 
ing this problem on the gravity side requires solving the 
full non-linear Einstein's equations, which typically can 
only be done numerically. Examples of such numerical 
studies include [1H5]. 

One of our purposes is to show that the problem on the 
gravity side can be simplified, at least in certain circum- 
stances. Inspiration comes from the so-called 'close limit 
approximation' (CLA) in the context of black hole 
mergers in four-dimensional general relativity in asymp- 
totically flat spacetime. The CLA is the statement that, 
once a single horizon forms around the two incident black 
holes, its subsequent evolution is well described by the 
Einstein's equations linearized around the final equilib- 
rium black hole, despite the fact that the initial horizon 
may or may not seem to be a small perturbation of the 
final one. In particular, the form of the gravitational ra- 
diation emitted to infinity in the merger-plus-ring-down 
phase is well described by the CLA llj. 

Following [5], we study isotropization of a homoge- 
neous plasma in a four-dimensional conformal field the- 
ory (CFT) in flat Minkowski space; on the gravity side 
this means that we work in the Poincare patch of AdSs. 
Note that, because of the homogeneity, the isotropization 
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FIG. 1; (a) Setup of Ref. [5]. (b) Our setup. The initial state 
is specified on a t = const, surface spanned by the radial 
coordinate r — see eqn. 



process involves exclusively non-hydrodynamics modes. 
Ref. jS] 'creates' a far-from-equilibrium state by acting 
on the CFT vacuum with an external anisotropic source 
— see Fig. [T] In contrast, we study the isotropization 
of a large number of anisotropic initial states in the ab- 
sence of external sources (see [7| [T^] for related work) . 
Each state is specified on the gravity side by an entire 
function on an initial-time slice (ITS), and hence it is 
characterized by an arbitrary number of scales. 

In this Letter we focus on the time evolution of the ex- 
pectation value of the gauge theory stress tensor; other 
observables will be considered elsewhere. Conservation 
of the stress tensor for a homogeneous plasma in the ab- 
sence of external sources implies that the energy den- 
sity £ (but not the entropy density) must be constant. 
Consequently, the final equilibrium state is known with- 
out solving for the dynamical evolution. On the CFT 
side it is a homogeneous, isotropic plasma with temper- 
ature oc f^/'* and pressure On the gravity side it 
is a static, isotropic black brane with the same tempera- 
ture. This makes the linearization of Einstein's equations 
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around the final state particularly simple. The dynamical 
evolution shows that an event horizon (but not necessar- 
ily an apparent horizon) is already present on the ITS 
for each of the states we consider. By comparing on the 
gravity side the full numerical evolution with the linear 
approximation (LA), we will see that the latter predicts 
with surprising accuracy the time evolution of the CFT 
stress tensor (see [9] for related observations), in analogy 
with the prediction of the gravitational radiation at in- 
finity by the CLA. We emphasize that the applicability 
of the LA is not guaranteed a priori, since in general our 
initial states are not near-equilibrium states. 
2. Holographic modeL For a conformal SU{N^) 
gauge theory, conservation and tracelessness of the stress 
tensor, together with homogeneity and rotational invari- 
ance in one plane (assumed for simplicity) imply that the 
stress tensor can be written as 

r 1 

(T;.) - ^diag^f, V,it), V^it), V^{t)\ , (1) 
V,{t) ^^£- ^AV{t) , VAt) ^l£ + l^Vit) , (2) 

in terms of a single function AV = 'Pt — 'Pl that measures 
the degree of anisotropy. Accordingly, the dual metric 
can be written as ^ 



2dtdr - Adt^ + S^e 
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dx. 



T?e^d^l , (3) 



where A, E and B are all functions of time t and of the 
radial coordinate r. In the absence of CFT sources these 
are subject to the following boundary conditions near the 
AdSs boundary at r cxi: 
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As usual, the normalizable modes and 64 (i) are not 
determined by the boundary conditions but must be read 
off from a full bulk solution that is regular in the interior. 
These modes are dual to the expectation value of the 
stress tensor. For the specific case of SU{N^) N = A 
super Yang-Mills this relation is 



£ = -304/4 and AV{t) = 3b4,{t) 



(5) 



Although £ is constant in time, a physical temperature 
can only be assigned to the system once (near) equilib- 
rium is reached. In this regime £ = 37r'*r*/4. 

In the generalized Eddington-Finkelstein coordinates 
([3]) Einstein's equations take the nested form 

= E (il)' + 2E' i; - 2E^ (6a) 

= E(S)' + |(E'B + S'E) , (6b) 

= A" +3S'B - 12E'E/E2 -f 4, (6c) 

= E+ ^(B^i]- A'E) , (6d) 

= E" + ^B'^ E , (6e) 



where h' = drh and h = dth , ^ 
along ingoing and outgoing null geodesies 



]^Adrh are derivatives 
respectively. 



Eqs. (6a I- (6c) are dynamical equations, whereas eqs. (6d) 
and (6el are constraints. If the former hold, then (6d) 



and (6e) are satisfied everywhere provided they are sat- 
isfied near the AdS boundary and on the ITS, respec- 
tively. As is clear from the causal structure in Fig.jlj the 
dynamical equations together with the constraints deter- 
mine the solution in the region labeled 'dynamics'. We 
find this solution by numerically evolving the full Ein- 
stein's equations following the procedure outlined in [S]. 

is a constraint on the possible initial states 



Eqn. (6e 



because it relates two of the metric functions on the ITS. 
We choose B as the independent variable because it is 
directly related to the CFT anisotropy through eqn. ([5|. 
Thus each initial state is specified by a constant 04 and a 
function of the radial cordinate B{t — 0, r). Note that for 
positive E the constraint (6c) implies E" < 0, which in 
combination with the asymptotic behavior E ~ r means 
that E will vanish at some r > on the ITS. Preliminary 
explorations indicate that this is a curvature singularity, 
and we will come back to this issue elsewhere. In any 
case, for all the initial states which our numerical code 
was able to evolve in a stable manner, the region where 
E = was hidden behind an event horizon and hence it 
had no effect on the physics. 

The static black brane solution of ([6| dual to a plasma 
in perfect equilibrium takes the form 



A^r^l-rty), 



B = 0, (7) 



with the horizon located at r,, = ttT. Considering small 
fluctuations around these equilibrium values one finds 
that A and E are unmodified at linear order, whereas 
the _B-fluctuation obeys eqn. (6b I with E and A as in 
0. 

3. Results. We considered around 1000 initial states, 
most of them generated by a random procedure (to be 
explained elsewhere). For some of these an apparent 
horizon was present on the ITS, for some others it was 
not. On the one hand we determined the time evolu- 
tion of each state by means of the full, non-linear Ein- 
stein's equations. On the other hand we solved the linear 
equation for B. In each case we read off the pressure 
anisotropy by extracting bi{t) from the near-boundary 
behavior (|4| . Fig. [2^ shows the result from the full Ein- 
stein's equations for a representative initial state. Fig.[2]3 
shows the difference between the full solution and the LA 
for this state. The ratio in the overall scales of the plots, 
2/10, gives a rough estimate of the accuracy of the LA, 
namely 20%, which is remarkable given that the evolution 
is definitely far-from-equilibrium. This follows from the 
thick blue curve in Fig. [2^, which shows that the pressure 
anisotropy is almost an order of magnitude larger than 
the energy density at some points during the evolution. 

Although there is no precise definition of entropy den- 
sity far from equilibrium, it is interesting to examine the 
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FIG. 2: (a) Solution B(t,z) (with z = 1/r) obtained from the full Einstein's equations. The initial profile B{t = 0,z) = 
^{z/zh)^ sm{8z/zit) is shown as a thick red curve. The thick blue curve shows AT'{t)/£ as obtained from the full Einstein's 
equations. The thin magenta curve shows the value of A'P(t)/£ as obtained from the linear approximation, (b) Difference 
between the full solution and the linear approximation. 
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FIG. 3: Time evolution of the areas of the event (top, blue) 
and apparent (bottom, red) horizons for the initial state of 
Fig. [2^. The red dot at the origin signifies that there is no 
apparent horizon for this state at the initial time. From that 
time until the start of the red curve there is no apparent 
horizon within the range of the radial coordinate covered by 
our grid, but there could be one at a deeper position. 



time evolution of the area densities of the event and ap- 
parent horizons, since these coincide with the entropy 
density in equilibrium. Fig. [3] shows that both of these 
quantities are larger at the end of the evolution than at 
the beginning, suggesting that entropy is indeed gener- 
ated. Incidentally, note that no entropy is produced in 
the LA, since A and S are unmodified. 

We define the isotropization time t^^^ as the time be- 
yond which AV{t)/£ < 0.1. Fig. |4] shows the isotropiza- 
tion times obtained from the full evolution of the 1000 
initial profiles that we considered, as well as the difi^cr- 
ences between the values of these times as determined by 
the full Einstein's equations and by the LA. We see that 
the LA works with a 20% or better accuracy for most 
states, and also that isotropization times are t^^^ < l/T, 
with T the final temperature. 

4. Discussion. Small perturbations around an equili- 
brated plasma can be described in linear response theory. 
Equivalently, perturbations around the dual horizon can 
be described by linearizing Einstein's equations around 




FIG. 4: Results for the isotropization times obtained from 
the full evolution of 1000 initial states, and for the difi'erences 
between the full and the linearized evolution (normalized by 
the full isotropization time). The height of each bar indicates 
the number of states in each bin. 

the equilibrium black hole solution. For a homogeneous 
but anisotropic perturbation one may expect the LA to 
be applicable whenever AV/£ ^1. It is remarkable 
that, in a strongly coupled CFT with a gravity dual, the 
LA actually works relatively accurately well beyond this 
limit. 

In Fourier space, one may distinguish between hydro- 
dynamic modes (HDMs), with dispersion relations uj{q) 
such that a; — ^ as 5 — 0, and quasinormal modes 
(QNMs), for which uj{0) ^ 0. If the perturbation is 
anisotropic but homogeneous then the relaxation back 
to equilibrium involves exclusively the QNMs. In this 
sense the dynamics we have studied can be thought of as 
the far-from-equilibrium dynamics of the QNMs. If the 
perturbation is small then these modes evolve towards 
equilibrium linearly, independently of each other and on 
a time scale set by the imaginary parts of their frequen- 
cies. One could imagine extending the description to not- 
so-small perturbations by including non-linearities in the 
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form of interactions between the QNMs, but naively one 
would expect this expansion to break down for order- 
one anisotropics. Instead, our results imply that, for 
a homogeneous, strongly coupled plasma with a gravity 
dual, the isotropization process is still reasonably well 
described by QNMs that evolve approximately linearly 
and independently of each other, even in the presence of 
large anisotropics (see [9j for related observations). This 
means that, just as in the near-equilibrium case, the re- 
laxation towards equilibrium is characterized by a few 
quasinormal frequencies. In fact, we have verified that 
the time evolution of the stress tensor is well described 
by expanding and evolving B in terms of a sufficient num- 
ber of QNMs. In particular, a naive (under) estimate of 
the isotropization time can be obtained from the imagi- 
nary part of the lowest quasinormal frequency, which is 
Imwo ~ 8.5T. Our initial states have anisotropics of or- 
der 1 < AV/S < 20, which gives 0.27 < Tt,,, < 0.62. 
The reason why this is an underestimate is that the de- 
gree of anisotropy carried by each individual QNM can be 
much larger, typically as large as AV/S ^ 500, the total 
anisotropy being much smaller due to cancellations be- 
tween different modes. With this value one gets Tt^^^ ~ 1. 

It would be interesting to understand the reason be- 
hind the relative accuracy of the LA in our setup. This 
will certainly involve comparing the next-to-linear term 
in the expansion with the linear term. Yet, demanding 
that the former is much smaller than the latter might 
yield too restrictive a condition, since subsequent terms 
in the expansion might partially cancel each other. 

We stress that the conclusions in this Letter hold for 
the purpose of computing the time evolution of the expec- 
tation value of the CFT stress tensor. Other observables, 
e.g. those considered in HSj, may or may not be well 
described by the LA. This is under investigation. 

Another interesting question is whether the LA may be 
applicable in less symmetric situations, for example in the 
absence of homogeneity, in which case HDMs will play 
an important role. Cases of particular interest include 
expanding plasmas, a simple example of which is a boost- 
invariant plasma. In this case the late-time solution is 
an expanding fluid with temperature (at leading order in 
the hydrodynamic expansion) T{t) = A/(At)-'^/^, with 
T the proper time and A a constant. The dual gravity 
solution is a black brane with a time-dependent horizon 
located at rh(T) = ttT^t) [14]. Thus it may be useful 
to linearize the Einstein's equations around the late-time 
solution (possibly including known gradient corrections) 
if A can be determined in terms of initial data defined at 
r > 0. This is currently under investigation. 

We close by stressing that the accuracy of the LA is 
a statement about the dynamics of black hole horizons. 
The LA is not expected to be applicable to the descrip- 
tion of strong gravitational dynamics in the absence of 
horizons. In particular, it is not expected to be able to 
describe the formation of a horizon — but it might be 



very useful in order to describe its subsequent evolution. 
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